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Abstract. We consider driven systems where the driving induces jumps in energy space: 
(1) particles pulsed by a step potential; (2) particles in a box with a moving wall; (3) particles 
in a ring driven by an electro-motive-force. In all these cases the route towards quantum-classical 
correspondence is highly non-trivial. Some insight is gained by observing that the dynamics in 
energy space, where n is the level index, is essentially the same as that of Bloch electrons in a 
tight binding model, where n is the site index. The mean level spacing is like a constant electric 
field and the driving induces long range hopping oc l/(n — m) . 



1. Introduction 

Consider a system which is described by a Hamiltonian H(X(t)), where the parameter 
X(t) is time dependent. For such system the energy E is not a constant of the 
motion. Rather, the driving induces spreading in energy space. Assuming that the 
system is prepared at t = in a microcanonical state, one wonders how the energy 
distribution pt(E) looks like at a later time. In particular, one may wonder whether the 
quantum pt{E) is similar to the corresponding classical distribution. In the "quantum 
chaos" literature it is customary to distinguish between a classical time scale r c / and a 
quantum breaktime t*. The latter goes to infinity in the u h — > 0" limit. A prototype 
model is the "quantum kicked rotator" pQ where the energy spreading is diffusive up 
to t* while for larger times one observes saturation due to a dynamical localization 
effect. 

In this work we analyze much simpler systems where the breaktime t*, if exists, 
is much larger than any physically relevant time scale. In fact one may assume that 
the time t of the evolution is comparable with the classical (short) time scale. In such 
circumstances one naively would expect quantum to classical correspondence (QCC). 
But in fact the theory is much more complicated 12] ■ One has to distinguish between 

• Detailed QCC 

• Restricted QCC 

Detailed QCC means that all the moments r — 1,2, 3,... of the quantum mechanical 
distribution pt(E) are similar to the classical result, while restricted QCC refers only 
to the r = 1, 2 moments. It turns out that the latter are very robust, while the higher 
moments (r > 2) might be much larger in the quantum case. Our first challenge would 
be to find and to analyze the worst case for QCC, for which all the r > 2 moments are 
classically finite but quantum mechanically divergent. We would like to see whether 
in such circumstances restricted QCC for r = 1,2 survives. 
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For completeness of this Introduction we summarize in Appendix A the reason 
for the robustness of restricted QCC. Our interest in QCC is motivated by the wish to 
develop a better understanding of driven systems. We would like to explore examples 
where QCC is far from obvious even for short times. In what follows we address 
4 problems that in first sight look unrelated: 

(1) Particles that are pulsed by a step potential (Fig.l) 

(2) Particles in a box with a moving wall (Fig. 2) 

(3) Particles in an electro-motive-force (EMF) driven ring (Fig. 3) 

(4) Wavepacket dynamics of Bloch electrons in a constant electric field 

In fact we are going to see that problems (l)-(3) share a common feature: In the 
classical description the energy absorption is associated with abrupt jumps in phase 
space. These jumps arc reflected in the quantum dynamics as a strong diffraction 
effect. This diffraction, which takes place in energy space, is the worst case for Bohr's 
QCC. It turns our that problem (1) can be solved exactly, while problems (2) and (3) 
reduce essentially to problem (4). Namely, the dynamics in energy space, where n is 
the level index, is essentially the same as that of Bloch electrons in a tight binding 
model, where n is the site index. The mean level spacing is like a constant electric 
field and the driving induces long range hopping oc l/(n — m) . This tight binding 
problem has an exact solution. The objectives of the present work are 

• To highlight the route towards QCC in the case of diffractive energy spreading. 

• To provide solutions and numerical demonstrations to the prototype problems. 

• To shed new light of the EMF-driven ring problem. 

• To illuminate the limitations of linear response theory in the mesoscopic context. 

The paper is structured accordingly. 

A few words are in order regarding the literature. The quantum treatment of 
the "moving wall" problem has started with Refs.|3J^, that were aimed in finding 
the steady state solutions for an expanding well. The interest in this model has 
further evolved within the study of the Fermi acceleration problem where the 
wall is oscillating. Recently the non-trivial features of the parametric [Hj and of the 
time-dependent wavepacket dynamics were illuminated. In the latter publication 
a satisfactory mathematical treatment of the non-stationary dynamics has not been 
introduced. Also the problem of Bloch electrons in a constant electric field has a long 
history. The concept of a Stark ladder was introduced by Wannier [H] to describe 
the energy spectrum of a periodic system in an electric field. Since that time it has 
become the subject of controversy P CHI IIH H21 E31 IHI CHI ESI EZI CHI CHI Eventually 
it has been realized that the electric field localizes the motion of the electrons, and 
induces a periodic oscillatory motion. 

2. Energy jumps in phase space 

If a Gaussian wavepacket is moving in a smooth potential, then its Wigner function 
evolves in a smooth manner which favors detailed QCC. But we would like to consider 
the "worst case" for QCC. Let us assume that the particle is prepared with some 
initial momentum p. This means in practice a very extended wavepacket with a very 
small dispersion in momentum. We turn on at t — a step of height V Btep . After a 
short time t we observe that the classical phase space distribution is torn into three 
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pieces (see Fig.l): phase space points that remain on the left side of the step; phase 
space points that have crossed the step from left to right; and phase space points that 
were all the time in the right side of the step. The jump in the kinetic energy of those 
points that have crossed the step is 

se c1 = -y stop (i) 

Classically we have in phase space points that move with the original kinetic energy 
and another set of points that have gone through an abrupt change of kinetic energy. 
Thus the energy distribution consists of two delta peaks. We would like to know what 
is the corresponding energy distribution in the quantum mechanical case. 

A similar phase space picture emerges in the analysis of the "moving wall" 
problem. As illustrated in Fig. 2 we have a particle of mass m and energy E 
bouncing back and forth inside a one dimensional box. One wall of the box is 
displaced with a velocity V^ all , which is assumed to be much smaller compared with the 
velocity v E — (2E/m) 1/>2 of the bouncing particle. Consider an initial microcanonical 
distribution. After a short time t some of the phase space points collide with the 
wall which is moving with velocity V mn - Consequently their velocity undergoes a 
change v i— > — v + 2V^ all , and accordingly the energy jump is: 

SE cl = -2mwBK.il (2) 

Thus after a short time the energy distribution consists of two delta peaks: one 
corresponds to those phase space points that did not collide with the moving wall, and 
the other corresponds to those phase space points that did collide with the moving 
wall. We ask what is the corresponding quantum result. Namely, how the probability 
is distributed among the energy levels in the quantum mechanical case. It is implicit 
that we are going to work in the adiabatic (wall location dependent) basis, else the 
question is mathematically ill defined. 

Possibly the most interesting and experimentally relevant model is that of a one- 
dimensional EMF-driven ring (Fig. 3). The classical analysis for this problem is very 
simple: each time that the particle crosses the EMF step its energy changes by 

SE cl = eV EMF (3) 

So also here we have energy jumps. Surprisingly this problem is interesting even if we 
do not add a scatterer. 

3. Beyond the Fermi golden rule, the semiclassical regime 

Both in the case of the "moving wall" and in the case of the driven ring we have after 
a short time a finite probability to find the system with a different energy. So we may 
say that there is some finite probability to make a transition 

E i — ► E + SE cl (4) 

Going to the quantum mechanical problem we may wonder whether or how we get 
from the Schrodingcr equation such transitions. We are used to the Fermi golden rule 
picture of transitions 

E i — ► E + 'W (5) 

where u> is the frequency of the driving. But here we do not have periodic ("AC") 
driving but rather linear ("DC") driving. Moreover, 5E cl is an h- independent quantity. 
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It turns out that indeed there exists a regime where the dynamics is classical-like 
(Fig. 4). This semiclassical regime is defined by the obvious condition 

5E cl » A (6) 

where A is the level spacing. In the case of the "moving wall" problem this condition 
can be written as 

V^ n » (?) 
mL 

where L is the size of the box. It is easily verified that this condition is just the 
opposite of the adiabatic condition. The case of the EMF-driven ring is somewhat 
richer. The condition that defines the semiclassical regime becomes 

V EMF » ^ (8) 

where L is the length of the ring. It is easily verified that this condition is just the 
opposite of the diabaticity condition. The diabatic regime is defined as that where 
transitions between energy levels of a "free" ring can be neglected. If there is a small 
scatterer inside the ring a stronger condition than diabaticity is required in order to 
maintain adiabaticity: 

VW«(l-fl)^ (9) 

where g ~ 1 is the transmission of the scatterer. The adiabatic regime is defined as that 
where transitions between the actual energy levels of the ring can be neglected. This 
is the regime where the Landau-Zener mechanism of transitions at avoided crossings 
(201 E] becomes significant. The three regimes in the EMF-driven ring problem are 
illustrated in the diagram of Fig. 5. 

Our main interest is in the non-trivial semiclassical regime as defined by Eq.©. 
In order to reconcile our semiclassical intuition with the quantum Fermi Golden rule 
picture we have to assume that the quantum dynamics self-generates a frequency 
"fa*;" = SE cl . Indeed it has been argued in Ref.0 that the non-perturbative mixing 
of levels on the small energy scales generate this frequency, while the re-normalized 
transitions on the large (coarse grained) energy scales are FGR-like. However, an 
actual mathematical analysis of the dynamics has not been introduced, and was left 
as an open problem. 



4. Particle pulsed by a step 



The simplest example for a semiclassical energy jump is provided by the "step 
problem" . The time dependent Hamiltonian is: 

h = 4— + i v. ; (10) 



r i o 

2m + \ V stap 

For this Hamiltonian "energy space" is in fact "momentum space", so it is more 
natural to refer to "momentum jumps" . Obviously we can translate any small change 
in momentum to energy units via SE — v E Sp, where v E = (2E/m) 1 ^ 2 is the velocity of 
the particle in the energy range of interest. 

The phase space dynamics after kicking an initial momentum state po at t = is 
illustrated in Fig. lb. It is clear that the emerging momentum distribution is 

'"' ■" Po)+^6[p-(po + 5 Pcl )) (11) 



Ptip) 



1 



Sip 
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where Sp d = —V stop /v E , and L is the spatial extent of the wavepacket. From here 
on we set L — 1 as implied by the standard density normalization of the momentum 
state e lp ° x . It is implicit in the following analysis that we assume a very extended 
wavepacket (v E t <C L). The emerging momentum distribution can be characterized 
by its moments with respect to p = p - Namely: 

((p-PoY) = Sp r ol xv E t r = 1,2,3,4,... (12) 

All the moments are finite and grow linearly with time. Below we are going to derive 
the quantum result. Omitting the trivial S(p — po) term, and the back reflection term, 
the final result for the forward scattering is 



fin 2 

Pt (p) = \(p\U\ Po )\ 2 = {p p Q)2 vlt 2 sine 2 



\{P~ (Po + Sp cl ))v E t 



(13) 



for which 

{Sp cl x v E t — sm(5p cl v E t) for r = 1 

Spl x v E t for r = 2 (14) 

oo for r > 2 

Let us compare the energy distribution in the classical and quantum-mechanical cases 
(Fig. 7). As the time t becomes much larger than h/V atop the semiclassical peak is 
resolved. But we never get detailed QCC, because all the high (r > 2) moments of 
the distribution diverge. 

It should be appreciated that the power law tails that we get here for the energy 
distribution are the "worst case" that can be expected. They emerge because the 
phase space distribution is torn in the momentum direction. In space representation 
this reflects a discontinuity in the derivative of the wavefunction. This explains why 
the tails go like 1/p 4 . We are going to encounter the same type of power law tails also 
in the other examples. 

4-1. Derivation of the quantum result: 

The rest of this section is devoted for the derivation of the quantum result and can be 
skipped in first reading. The momentum states are denoted as \p). In order to simplify 
the calculation we approximate the dispersion relation, within the energy window of 
interest, as linear E = v E k. This implies that back-reflection is neglected. Once the 
step is turned on, \p) are no longer the stationary states. The new stationary states 
are 

\k) .— » &{-x)c lkx + Q{x)S k+u ^ x (15) 

where we use the notation u = Sp cl . Note that these form a complete orthonormal set 
in the sense (fci|fc 2 ) = 27r<5(fci — fc 2 ). The transformation matrix from the old to the 
new basis is 

e- i(p - k)x dx+ / e - l(p - k - u)x dx (16) 

-oo JO 

% % 

= nS(p - k) H + n5(p - k - u) (17) 

p—k p—k—u 

Before we go on with the calculation we note that the following elementary integral 
can be found in any mathematical handbook: 

r°° dk i 



g M = T{ -(e^*-e Vlt ) pifri, t>0 



2-rr (k - p 2 )(k - pi) 2(p 2 -pi) 
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We notice that the result on the RHS if finite for p 2 = p\ while in fact it should 
diverge. This suggests that there is a missing delta term Ce lp2t 5(p2 — Pi) where C is 
a constant. In order to find this constant we have regularized this Fourier integral: 



dk k- p 2 k- pi _ iM 



/oo 
e 
2 7 r(fc-p 2 ) 2 + ( 5 2 {k- Pl ) 2 +S 2 



(e^t _ e ipit) + J e ip2t (5(p 2 - Pl ) t>0 



2( P2 -Pi) y ' 2 

With the above we can calculate the matrix elements of the evolution operator: 

(p\U\ Po ) = V e-* B ** (p\k) (k\ Po ) = — / e~ iv Ekt (p\k) (k\ Po ) dk 
= tt5{ P - p ) (e lVEpt + e lVE( -P +u ^ 

(p -Po)(p -Po - u)\ ' 
We have (p\U(t = Q)\po) = 2nS(p — po) as required. The interesting part of this 
expression is the second terms which is non vanishing for p ^ pq. Taking its absolute 
value and squaring we get after some algebra Eq. . 

5. Particle in a box with a moving wall 

5.1. The Schrodinger equation: 

We consider a particle in an infinite well. The left wall is assumed to be fixed at 
x = xq, while the right wall at x = X(t) is moving with constant velocity X = V mn - 
The size of the box is L(t) = X(t) — xq. Classically the dynamics is very simple: each 
time that the particle hits the moving wall its energy jumps by 5E cl = 2mvV„ s . n . In 
the quantum mechanical case we work in the adiabatic basis. The adiabatic energy 
levels and the eigenstates for a given value of L are: 

1 firh \ 2 

E n = w- -r n ( 19 ) 



2m V L 

¥ n \x) = (-ly j - sin f — x) (20) 



L V L 

We use the standard prescription in order to write the Schrodinger equation in the 
adiabatic basis. Using the notations of Ref.|22] the equation is written as 

—jr = -jE n a n + iX 2 t A nm a m (21) 
m 

where 

A nm = i(* (n) \^ (m) ) (22) 

Hence, the Schrodinger equation for the problem in the adiabatic basis is O IB] : 
da n i K-aii 2nm 

— = -- En a n -— ^ -j— a m (23) 
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5.2. The generated dynamics: 

Let us assume that the initial preparation is a„(0) = S nm . The mean level spacing 
for the ID box is A = irhv E /L. If 5E cl <C A one finds out, by inspection of Eci. (|23|) . 
that the dynamics is adiabatic, meaning that a n (t) ~ 5 nm . On the other hand, if 
5E cl A, one expects to find a semiclassical transition E i— > E + SE cl . 

How can we explain the E i— > E + 5E& transition from quantum-mechanical point 
of view? For this purpose we can adopt the core-tail picture of Ref. [231 : The 'core' 
consists of the levels that are mixed non-perturbatively; The 'tail' is formed by first 
order transitions from the core. Originally this picture has been applied to analyze the 
energy spreading in "quantum chaos" driven systems. Here, the (non-chaotic) moving 
wall problem allows a much simpler application \7\ . The analysis is carried out in two 
steps which are summarized below. 

The first step in the "core-tail" picture is to analyze the parametric evolution 
which is associated with Eq. (J23J) ■ This means to solve Ea. (|23[l without the first term 
in the RHS. (This is the so-called sudden limit). Obviously the resultant a n (t) is a 
function of 5X = V wan i, while K™n by itself makes no difference. The solution depends 
only on the endpoints x(0) and x{t). By careful inspection of Eq.J2HJ) one observes 
that a level is mixed with the nearby level whenever the wall is displaced a distance 
A E /2, where A E = 2irh/(mv E ) is the de Broglie wavelength. The time scale which is 
associated with this effect is obviously 

= ±£ (24) 

''wall 

The second step is to analyze the actual time evolution. This means to take into 
account the effect of the first term in the RHS of Eq. l|23|) . and to understand how 
the resultant a n (t) differs from a n (t). One observes that the 'parametric' mixing of 
nearby levels modulates the transition amplitude. The modulation frequency is 

V - — (25) 

T qm 

Once combined with the FGR Ea. (|25(l it leads to the anticipated semiclassical result 
Eq.©. It is not difficult to argue that the period of this semiclassical transition is 

r, = - (26) 

We 

which is the time to make one round between the walls of the well. Since we are 
dealing with a simple ID system this coincides with the Heisenberg time: 

tu = = r cl (27) 

The ratio T cl /r qm determines the number of nearby level transitions per period. 
Obviously the semiclassical condition Eq.J7J) requires this ratio to be much larger 
than unity. The disadvantage of the above heuristic picture is that it does not lead to 
a satisfactory quantitative results. Therefore, in later sections we discuss an optional 
route of analysis via a reduction to a tight binding model. 

5.3. Numerical Simulation 

The solution of Eq.JS^J) becomes very simple if we make the approximation 
L{t) w Lq. This holds as long as the wall displacement is small. We have 
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verified that the associated numerical error is very small. Using units such 
that L = m = h = 1 we define a diagonal matrix E = diag{7r 2 n 2 /2} and a non- 
diagonal matrix, W — {~i2anm/ (n 2 —m 2 )} with zeros along the diagonal, and where 
a = V mU /L. The evolution matrix in the adiabatic basis is obtained by exponentiation: 

U(t) = cxp [-i t (E + W)} (28) 

Fig. 8 illustrates the time dependence of probability distribution |a„(i)| 2 for a particle 
initially prepared at no — 50. Fig. 8a displays the solution in the adiabatic regime: 
the particle stays at the same level. Fig. 8b displays the solution in the semiclassical 
regime: at each moment the particle partially stays at the same energy, and partially 
makes classical-like transition to the next energy strip. Fig. 9 highlights the energy 
splitting of the wavepacket during the transition. 

If we want to avoid the L(t) Lq approximation, the price is a time dependent 
E and W matrices. Then the calculation should be done in small dt time steps: 
( 

U(t) = Y[ exp [-i dt W(t')] exp [-t dt E(t')} (29) 

t'=dt 

The state of the system is described by a truncated column vector a = {a n } of 
length N. Optionally it is possible to represent the state of the system in the Fourier 
transformed basis. The elements of the Fourier transformed vector are labeled 
by k — (2Tt/N)h, where n mod(-/V) is an integer. The practical implementation of 
Ea. H29|) is greatly simplified if W nm is a function of the difference n — m. In such case 
W is transformed into a diagonal matrix W. Consequently one can use the standard 
fast Fourier transform (FFT) algorithm in order to propagate a given state vector. 
Namely, 

t 

a(t) = Yl FFT" 1 exp [ - i dt W{t')] FFT exp [ - i dt E(t')] a(0) (30) 

t'=dt 

where both E and W are diagonal. In the moving wall problem W nm is mainly 
proportional to l/(n— to), so the FFT method is applicable if we restrict the energy 
range of interest. In the next section we shall consider the EMF-driven ring problem, 
leading to a very similar evolution equation, where the FFT method is strictly 
applicable. 



6. Particle in an EMF-driven ring 

6.1. The Schrodinger equation: 

We consider a ID-ring driven by an EMF (Fig. 3). The EMF is induced by a time- 
dependent flux which is described by the vector potential 

A(x,t) =®(t)5(x-x ) (31) 
This means that the electric field is 

£{x) = V EMF 6{x - x ) (32) 
where V EMF = — $ = const. The Hamiltonian that generates the dynamics is 

Umt)) = ^{p--A{x,t)) 2 (33) 
with periodic boundary conditions over x. The length of the ring is L. 
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Classically the dynamics is very simple: each time that the particle crosses x — Xq 
its energy jumps by 5E cl = eVE MP . In the quantum mechanical case it is convenient 
to work in the so-called diabatic basis. The diabatic energy levels for a given value of 
$ are 

where $o — ^irhc/e. See Fig. 4. The Schrodinger equation that describes the time 
evolution in the diabatic basis is found using the same procedure as in the case of the 
moving wall. We have to find the A nm as defined in Ea. (|22l) where now X = <&. The 
only extra difficulty is in finding the eigenstates of Ea. (|33|l because A(x) depends 
on x. The calculation becomes much simpler if we realize that they are related by a 
gauge transformation to the eigenstates of a much simpler Hamiltonian: 

Namely, 

¥ n \x) = exp(^A(x)j¥ n \x) (36) 
/ ie , ,\ 1 / 2im \ 

= exp I he ^ ) x "Tl exp \~ x ) ^ ' 

exp [v^- (—+n)x) (38) 



where in the last line we set Xq = and the gauge function is 

A(x) = y x (39) 

Using the above result we get 
i 1 

A nm = -- (40) 

9q n — m 

and accordingly 

dd n i ^EMF \ ^ 1 / a t \ 

^r = -H Enan + ^r ^ ^ flm (41) 

6.2. The generated dynamics: 

The dynamics of an EMF-driven ring is very similar to the dynamics in the moving 
wall problem. This is obvious from the phase space picture, and also by inspection 
of the equation for a n (t). Also the "core tail" heuristic picture of section 5.2 is easily 
adapted. The parametric scale that signifies mixing of nearby levels is now 5X = $o 
instead SX = A E /2 leading to the quantum time scale 

Tqm = v — 

The classical period is 

r OJ = - (43) 
and the semiclassical condition can be written as r am <C t cI . 



(42) 
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Since the energies are time dependent we have to use Eq. I|29|l for the calculation 
of the time evolutions. Furthermore, W is diagonal in the momentum representation, 
and therefore we can use the FFT method Eq. with W = diag{— a (k — tt)}, 
where k = (2ir/N)n is defined mod(27r). The results of the simulations are presented 
in Fig. 10a and Fig. 11. We shall further discuss these results in the next sections. 



7. Bloch electrons in a constant electric field (I) 

If we focus our interest in small energy interval, then in both cases (moving wall, driven 
ring) the Schrodinger equation in the adiabatic basis is approximately the same as that 
of an electron in a tight binding model, where n is re-interpreted as the site index: 

^JT = ~iE n a n + a V" — - — a m (44) 
at * — ' n — to 

with E n = en. We use from here on H = 1 units. The scaled rate of the driving a 
is re-interpreted as the hopping amplitude between sites, while the levels spacing e is 
re-interpreted as an electric field. Assuming that the electron is initially at the site 
no we would like to find out what is the probability distribution 

Pt{n) = \a n {t)\ 2 (45) 

It is obvious that the adiabatic regime a <C e corresponds to a large electric field that 
localizes the electron at its original site. In the other extreme (a > e), if the effect 
of e could have been ignored, we would observe unbounded Bloch ballistic motion. 
The effect of finite e is to turn this motion into Bloch oscillations. We shall find below 
that the electron performs periodic motion which we illustrate in Fig. 6: While the 
wavepaket drifts with the electric field to the right, it shrinks and disappears, and at 
the same time re-emerges on the left. If we run the simulation as a movie, it looks as if 
the motion is from left to right. Still it is bounded in space due to this "re-injection" 
mechanism. 

First of all we solve the equation for e = 0. The Hamiltonian is diagonal in the 
momentum basis k and therefore the general solution is 

a n (t) = A k e ll - kn ~^ = f —A k e^ kn - Uk ^ (46) 
The dispersion relation is found by transforming the Hamiltonian to the k basis: 

„ — ik(n— m) 

LOk = ice y — a [tt — k] (47) 



-ikno 



If we place at t = an electron at site no, then a n = <5 n . ri0 , and hence Ak = c 
Then we get 

, . sin7rat , 

a n {t) = ~ T — r (48) 

ir{at + n — no) 

Turning to the general case with e ^ we substitute a n (t) = c n (t)e~ iEnt and get the 
equation 

J i(n-m)et 

=*■ = <* E c ™ ( 49 ) 

at * — 4 n — m 
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This more complicated equation is still diagonal in the k basis: 
dC k 



dt 



= -iu) k (t)C k 



where 

LOk{t) = a [it — mod(fc + et, 2ir)] 
and its solutions is 



(50) 
(51) 



Cfc(t) 



C fc (t=0) exp 



i / uo k {t')dt' 



o 



= < 



-ifcn +ic(((fc-7r)t+^|-J 

-ikn a +ia((k-K)t+^- + ^^j 
-ikn +ia((k+v)t+^-) 



for e > and < k < k 
for £ > and k < k < 2ir 
for £ < and < k < k 
for £ < and k < k < 2ir 



which is valid for < t < 2tt /\e\ and should be continued periodically in time. We have 
used the notation k = —et mod(27r). Now we can go back to position representation: 



c n {t) 



k ^dk 
2tt 



r~t Akn I 

G fc e + 



k(t] 



dk 
2tt 



ikn 



(52) 



Taking the absolute value and squaring we get the following result for the probability 
distributions: 



Pt 



2 sin 2 (l£t(n-n + a (t- ffi)) 
W ~ V £ ) (n-n + at) 2 (n - n + a(t - ^)) 2 



(53) 



The above formula is valid for < t < 2tt/\e\ and it should be continued periodically 
in time. Fig. 6 and Fig. 10a illustrate the dynamics both schematically and numerically. 
In the next section we further discuss the nature of this dynamics. 



8. Bloch electrons in a constant electric field (II) 

In order to appreciate the significance of the oc l/{n— m) hopping we solve again 
the problem of Bloch electrons in a constant electric field, but this time with the 
"conventional" nearest neighbor hopping: 



da n a 



with E n = en. The initial preparation at t—0 is a n — 5 n , r , 



-iE-t 



c n and get the equation: 



C n +1 



e l£t c„_!) 



dt 2 

This equation becomes diagonal in the k basis: 

—jj- = -iuj k (t)Ck 
where 

Wk(t) — asin(£t + fc) 



(54) 

We substitute 
(55) 

(56) 
(57) 
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Its solutions is 

rt 



C k {t) = C7 fc (t = 0) xexp 



i / uj k (t')dt' 



o 



(58) 



Solving the above integral and making the inverse Fourier transform we obtain: 



c n {t) = J n - na ( — sin ( \et ) ) (59) 



V s \2 

where J() is the Bessel function of the first kind. Taking the absolute value and 
squaring we get the probability distribution: 

2 



Pt(n) 



, 2a . (I 

Jn-n | — Sin I -St 



(60) 



Fig. 10b illustrates the dynamics. As in the previous problems we can distinguish 
between two time scales. One is related to the diagonal part of the Hamiltonian, and 
the other one to the hopping term. Keeping the same notations as in previous sections 
these are 

r cl - 2ir/e (61) 
r qm = 1/a (62) 

The nature of the dynamics in the case of oc 1/ (n—m) hopping and in the case of near 
neighbor hopping is quite different, as it can be appreciated by comparing Fig. 10a and 
Fig. 10b. This is related to the additional symmetries in the latter case. In order to 
explain this point let us use the notation U(a, s) that emphasizes that the evolution 
depends on two parameters, the first one is associated with the kinetic term W = w(p) 
and the other one with the potential term E — e(x). For clarity we use x for the 
position coordinate and p for the quasi-momentum. In both cases we have the anti- 
unitary symmetry (x,p) [x, —p), that maps E to E and W to —W. Consequently 
U(a;s) is mapped to U(a; —e). This implies that the spreading does not depend on 
the direction of the electric field. This is a peculiarity of tight binding models. The 
conventional time reversal symmetry, for which the kinetic term W is left invariant, 
is (x,p) i ► (x, ir—p). This symmetry characterizes the near neighbor hopping, but not 
the oc I /(n—m) hopping. This symmetry implies that the spreading looks the same if 
we reverse the signs of both a and e, which is like reversing the time. If we combine 
the two symmetries we deduce that the dynamics, in the case of the near neighbor 
hopping, should be indifferent to the sign of a. Note that the combined symmetry 
that leads to this conclusion is the unitary mapping (x,p) (x,p+n). Thus, in both 
cases [oc l/(n— to) hopping and near neighbor hopping] we have generalized Bloch 
oscillations, but in the former case they are unidirectional (Fig. 6), while in the latter 
case they are bi-directional. 



9. Discussion 



Within Linear response theory (LRT) the energy absorption of a quantum system is 
determined by the correlation function of the perturbation term. In general one can 
argue that there is a very good QCC for the correlation functions, and hence one 
expects restricted QCC in the energy absorption process. The persistence of restricted 
QCC in the t — > oo limit requires the additional assumption of having a coarse grained 
Markovian-likc behavior for long times. Depending of the context one should further 
assume that the environment supplies both weak decoherence effect that makes the 
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break time t* irrelevant, and a weak relaxation effect so as to achieve a steady state. 
Then it is possible to use the same argumentation as in the derivation of the central 
limit theorem in order to argue that all the higher moments become Gaussian-like. 

Thus, the common perception is that the leading result for the response of a driven 
system should be the same classically and quantum mechanically. For example, such 
is the case if one calculates the conductance of a diffusive ring [21] : The leading order 
result is just the Drude expressions, and on top there are weak localization corrections. 

The above reasoning illuminates that the long time response is based on the short 
time analysis. Moreover, one realizes that the second moment of the evolving energy 
distribution has a special significance. Still, all the above observations are within a 
very restrictive framework of assumptions. In practice it is of much interest to explore 
the limitations of LRT, and to obtain a more general theory for response. 

The theory for the response of closed isolated driven quantized chaotic mesoscopic 
systems is far from being trivial, even if the interactions between the particles are 
neglected. In the case of a generic quantized chaotic systems two energy scales are 
involved: the mean level spacing A oc h d , where d = 2,3 is the dimensionality of the 
system, and the semiclassical energy scale h/r c i. It is implied (see the mini-review 
of Ref. 25]) that there are generically three regimes depending on the rate X of the 
driving: 

• The adiabatic (Landau Zener) regime 

• The Fermi-golden-rule (FGR, Kubo) regime 

• The semiclassical (non-perturbative) regime 

Most of the literature in mesoscopic physics is dedicated to the study of the dynamics 
in either the adiabatic or the FGR regimes. The existence of a non-perturbative regime 
|2fil 12*5] is not yet fully acknowledged, though it has been established numerically in 
the RMT context [77|. 

Driven one-dimensional systems are non-generic because typically the semiclassi- 
cal energy scale coincides with the mean level spacing. In other words: the Heisenberg 
time tu = 2irh/A is the same as the classical time r c / rather than being much larger. 
Indeed we have seen that in the "moving wall" problem we have just two regimes: the 
adiabatic regime and the semiclassical regime. 

The EMF-driven ring is a prototype problem in mesoscopic physics. It is richer 
than the "moving wall" problem because a small scatterer introduces a very small 
energy scale, the level splitting, and hence we have three regimes rather than two: 
adiabatic, diabatic and semiclassical. 

The semiclassical regime in the study of EMF-driven rings has not been explored 
so far. One important observation is that contrary to LRT the gauge of the vector 
potential does matter. Most of past studies assume that the vector potential is 
A(x,t) = <&(t)/L. It is true that in LRT the same result for the conductance is 
obtained with A(x, t) — $(t)5{x — xq). If we try to go from A(x, t) to A(x, t) using a 
guage transformation, the "price" is a modified V(x) that features a linear ramp with 
a step- like drop at x = xq. This modification of V(x) can be neglected only in the 
LRT regime. The semiclassical condition of Eq.© is just that opposite of this LRT 
requirement. 

The semiclassical dynamics implies diffractive energy spreading. The mixing 
of levels in the small energy scales induces jumps in energy space. The realization 
that this diffractive energy spreading can be re-interpreted using a tight binding 
Bloch model follows in spirit the celebrated reduction pp of dynamical localization 
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in periodically kicked systems to a tight binding Anderson problem. An interesting 
feature is the hopping that goes like oc l/(n — to). This hopping leads to 
an unidirectional rather than bidirectional Bloch oscillations, as implied by the 
semiclassical reasoning. 



10. Appendix: The robustness of restricted QCC 

The simplest way to illuminate the robustness of the second moment is by adopting 
a heuristic phase space picture language. Given two operators A = A(x,p) and 
B = B(x,p), with Wigner Weyl representation A wv ,(x,p) and B ww (x,p) we have 
the exact identity 

~ ~ f dxdp 

trace(AB) = / A ww (x,p)B ww (x, p) (63) 

If we can justify the replacement of A ww (x,p) by A{x,p) and B ww (x,p) by B(x,p) 
then we get QCC. The rule of the thumb is that in order to justify such an 
approximation the phase space contours of A(x,p) and B{x 1 p) should be significantly 
different. Otherwise the transverse structure of the Wigner- Weyl functions should be 
taken into account. This reasoning can be regarded as a phase space version of the 
stationary phase approximation. 

Let A = [H(x,p)} r the rth power of the Hamiltonian H = H(x,p), and let 
B = p(Hq(x,p)) a stationary preparation with the Hamiltonian 7io = Hq(x,p). In 
such case trace(AB) is the rth moment (TV") of the energy. If H = Hq + XV, and A 
is not large enough, then we do not have detailed QCC. This is discussed throughly 
in Ref. 28 . But at the same time, irrespective of A, restricted QCC is robust. The 
reason is that for the first two moments we have the identities 

(H) = (H ) + (V) (64) 
(H 2 ) = (H 2 )+2(H )(V) + (V 2 ) (65) 

Thus the calculation of trace(AB) with A = [H(x,p)] r reduces to the calculation of 
trace(AB) with A = [V(x,p)] r . We assume that V and H are not related in any special 
way. It follows that we have robust QCC for all the moments of V, and consequently 
also for the first two moments of H, irrespective of A. 

In order to generalize the above reasoning to time dependent Hamiltonians, it is 
convenient to adopt the Heisenberg picture. Given that the system is prepared in a 
stationary state at t = 0, one can prove that 

(H(t) 2 ) - (H(0) 2 ) o = ((H(t) - n(0)) 2 ) (66) 

where TL{t) is the Hamiltonian TL{X{t)) in the Heisenberg picture. Such relation cannot 
not be generalized to higher moments because of lack of commutativity. Using 

f ^ f - **® «"> 

where V = dH/dX, we can express ((7~C(t) — Ti(0)) r )o as an integral over the 
correlation functions of the perturbation V(t). The QCC for these correlation 
functions is robust, and hence the QCC for the second moment is also robust. 
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Fig.l: (a) Left panel: Picture of the potential. Before t=0 the potential is zero (dashed line). 
After t=0 the potential is the step function (solid line), (b) Right panel: Phase space picture. Before 
t=0 there is no potential and the momentum is constant (dashed line). The piece of the distribution 
that has passed x=0 after t=0 is boosted with 5E cl = — Kte P . 
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Fig. 2: (a) Left panel: Potential well. The right wall is moving with a constant velocity y wa n. 
(b) Right panel: Phase space picture. If the wall were not moving, the distribution would evolve 
along the dashed line. If Vwaii is non zero, an energy jump 5E a] = — 2mv E V v , B ,u is associated with the 
collision, and one obtains the distribution which is illustrated by the solid line. 
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Fig. 3: (a) Left panel: Ring with EMF. (b) Right panel: Phase space picture. Without the EMF 
the momentum is a constant of the motion (dashed line). Else an energy jump 6E cl = eVsMF is 
associated with each crossing of the EMF step. The emerging phase space distribution is illustrated 
by the solid line. 
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L(t) 

Fig. 4: (a) Left panel: The energy levels of a one dimensional box as a function of its width L(t). 
For the purpose of comparison with other figures L(t) is decreasing as we go to the right (the box 
becomes smaller) so the levels are going up. The solid line illustrates adiabatic dynamics, while the 
"jumps" illustrate semiclassical dynamics, (b) Right panel: The energy levels of a one dimensional 
ring as a function of the Aharonov-Bohm flux. The solid line illustrates diabatic dynamics, while the 
"jumps" illustrate semiclassical dynamics. The dashed line illustrates adiabatic dynamics. 
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Fig. 5: The three regimes in the EMF driven ring problem. See text. 




Fig. 6: The unidirectional oscillations of Bloch electrons with oc l/(n — m) hopping: As the wavepaket 
slides to the right it shrinks, while being re-injected on the left, where it re-emerges. This should 
be contrasted with conventional bi-directional oscillations of Bloch electrons with nearest neighbor 
hopping. 
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Fig. 7: (a) Left panel: The classical energy distribution Eg, 1111 some time after a step potential 
is turned on. In this illustration V - Btcp <0. (b) Right panel: the corresponding quantum mechanical 
energy distribution calculated with Eg, 1131 . 
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Fig. 8: Density plots of the probability distribution as a function of time for the moving wall problem, 
(a) Left panel: Adiabatic regime. The probability stays at the same level. In order to clarify the 
connection with Fig. 4a we have added an E axis. The constant energy dashed lines are for guiding 
the eye. The populated adiabatic level goes up in energy which implies that the particle is steadily 
increasing its energy, (b) Right panel: semiclassical regime. The probability jumps in energy space. 
Note that with respect to the E axis we have the steps of Fig. 4a. The parameters of these simulations 
were L=m=7i=l and V wa u=0.l7r for (a) and V wa ii=57r for (b). Note that for n=50 the classical period 
is r c i=0.0127. The vertical dashed lines indicate two representative times t=T cl and t=1.5r c i. 
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Fig. 9: Plots of the probability distribution |a n | 2 at the two representative times that were indicated 
by the vertical dashed lines in the previous figure. In the classical limit the energy distribution 
consists of delta peaks instead of broadened peaks, in complete analogy with Fig. 7. 
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Fig. 10: (a) Density plots of the probability distribution as a function of time for the EMF driven ring 
problem. See the caption of Fig 8 for presentation details. The parameters are L = m = h = e = l 
with Vemf = 588840. Note that for n = 50 the classical period is t c i = 0.00021. This is 
approximately the same as solving the wavepacket dynamics for Bloch electron in electric field Eg. 1441 
with a = 93717 and e = 31416. (b) Wavepacket dynamics for Bloch electron in electric field with 
near-neighbor hopping. The parameters are m = h = e = 1 with a = 14139 and e = 3141.6. 
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Fig. 11: Plots of the probability distribution \a n \ 2 at representative times (as indicated). The solid 
lines are the solution of Eg. 1441 for Bloch electrons, while the dotted lines are the exact numerical 
solutions for the EMF driven ring, taking into account the quadratic (rather than linear) dependence 
of the eigen-energies on <J>. 



